Mon. Not. R. Astron. Soc. OOP. [THTSl (2010') Printed 4 July 2011 (MN 1^1^ style file v2.2) 

Galaxy And Mass Assembly (GAM A): the red fraction 
and radial distribution of satellite galaxies 

Matthew Prescott,^* I.K. Baldry/ P.A. James/ S.P. Bamford,^ J. Bland-Hawthorn,^ 
S. Brough,^ M.J.I. Brown, ^ E. Cameron,^ C.J. Consehce,^ S.M. Croom,^ 
't^'-S.F. Driver,^ C.S. Frenk,^ M. Gunawardhana,^ D.T. Hill,^ A.M. Hopkins,"^ 
(3 D.H. Jones, ^ L.S. Kelvin,'' K. Kuijken,^ J. Liske,^^ J. Loveday,^^ 
<N R.C. Nichol,i2 p Norberg,^^ H.R. Parkinson,!^ J.A. Peacock,!^ S. Phimpps,^^ 
3 K.A. Pimbblet,^ C.C. Popescu,!^ A.S.G. Robotham,^ R.G. Sharp,!^ W.J. Sutherland,!^ 
E.N. Taylor,^ R.J. Tuffs, ^^ E. van Kampen,!*^ and D. Wijesinghe^ 



O 



Oh- 



^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead, CH4I ILD, UK 
^ Centre for Astronomy and Particle Theory, University of Nottingham, University Park, Nottingham NG7 2RD, UK 
^Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia 
_ . , '^Australian Astronomical Observatory, PO Box 296, Epping, NSW 1710, Australia 
\^ ' ^School of Physics, Monash University, Clayton, Victoria 3800, Australia 

^Department of Physics, Swiss Federal Institute of Technology (ETH-Ziirich), 8093 Ziirich, Switzerland 
^School of Physics & Astronomy, University of St Andrews, North Haugh, St Andrews, KY16 9SS, UK 
I , ^ Institute for Computational Cosmology, Department of Physics, Durham University, South Road, Durham DHl 3LE, UK 
Q ■ ^Leiden University, P.O. Box 9500, 2300 RA Leiden, The Netherlands 
5— ( ' ^'^ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 
^y^ ' Astronomy Centre, University of Sussex, Falmer, Brighton BNl 9QH, UK 

C^ , ^^ Institute of Cosmology and Gravitation (ICG), University of Portsmouth, Dennis Sciama Building, Burnaby Road, POl 3FX, UK 
^''Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3H.I, UK 
' ^"^ Astrophysics Group, HH Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 ITL 
" ' ^^ Jeremiah Horrocks Institute, University of Central Lancashire, Preston PRl 2HE, UK 

^ , ^^ Research School of Astronomy & Astrophysics Mount Stromlo Observatory Cotter Road Weston Creek, ACT 2611 Australia 
' - ^"^ Astronomy Unit, Queen Mary University London, Mile End Rd, London El 4-NS, UK 

^^ Max Planck Institute for Nuclear Physics (MPIK), Saupfercheckweg 1, 69117 Heidelberg, Germany 



^ 



v^ ' Accepted 2011 June by MNRAS; received in original form 2011 January 

o ^ 



X 



ABSTRACT 

We investigate the properties of satellite galaxies that surround isolated hosts 

within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And 

H I Mass Assembly survey. Making use of isolation and satellite criteria that take into 

account stellar mass estimates, we find 3 514 isolated galaxies of which 1426 host a 
total of 2 998 satellites. Separating the red and blue populations of satellites and hosts, 
using colour-mass diagrams, we investigate the radial distribution of satellite galaxies 
and determine how the red fraction of satellites varies as a function of satellite mass, 
host mass and the projected distance from their host. Comparing the red fraction of 
satellites to a control sample of small neighbours at greater projected radii, we show 
that the increase in red fraction is primarily a function of host mass. The satellite red 
fraction is about 0.2 higher than the control sample for hosts with 11.0 < logj^o-^* < 
11.5, while the red fractions show no difference for hosts with 10.0 < logj^o ■^* < 10.5. 
For the satellites of more massive hosts the red fraction also increases as a function of 
decreasing projected distance. Our results suggest that the likely main mechanism for 
the quenching of star formation in satellites hosted by isolated galaxies is strangulation. 

Key words: surveys - galaxies: star formation - galaxies: interaction - galaxies: 
evolution - galaxies: formation - galaxies: dwarf. 
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1 INTRODUCTION 

In recent years satellite galaxies have received much atten- 
tion in both observational and theoretical studies in order to 
establish their role in the formation and evolution of galax- 
ies. In the current ACDM models of the Universe, galaxies 
are assembled in a hierarchical fashion, whereby small halos 
of dark matter (DM) merge to form larger halos, in which 
baryonic matter then cools and condenses to form stars. In 
this framework satellite galaxies are associated with sub- 
halos of dark matter residing within the virial radii of larger 
halos, which are believed to be left over from an earlier as- 
sembly phase of their host. Thus the measurement of the 
spatial distribution of satellites, both in terms of their an- 
gular and radial distributions, can provide an insight into 
the mass accretion histories of galaxies. 

The most common way to determine the radial 
distribution is to calculate the projected density of 
satellites surrounding samples of more luminous hosts, 
which requires redshifts and photometry of satellites and 
hosts. Previous studies which attempted to constrain the 
small - scale galaxy correla t ion function (ILake fc Tremaind 
198d; iPhillipps fc ShanQ Il987l : IVader fc Sandagd Il99ll : 



Lorrimer et al 



1994) or investigated the compan ions of field 



eUipticals IjMadore et "ai1l2004l : ISmith et al.ll2004 ) found that 
the projected density of satellites, as a function of radius, 
has a profile that can be described as a power law of the 
form E(_R) ex R" , with a slope, a ranging between —0.5 
and —1.25. Subsequently, more accurate measurements of 
the projected density have been made, making use of larger 
and more complete redshift surveys, s uch as the 2-degree 
Field Galaxy Redshift Survey f2dFGRS. [Colless et al.ll200ll ') 
and the Sloan Digital Sky Survey fSDSS. lYork et al.ll200(J ). 
allowing the production of large statistical samples of satel- 
lites and hosts selected using various we l l defined criteria 
(ISales & Lambas '2004'; 'Ch en et al.l l2006l : lAnn et al.l l2008l : 
[Chenii2008 ; BaiUn et al. 20oi v 

Studies which have examined the radial distribution of 
satellites as a function of host luminosity, colour or morphol- 
ogy, have obtained mixed results. The first to divide their 
sample into early and late types was lLorrimer et al.l ([199J), 
who found the projected density profile of the satellites of 
early-type hosts to have steeper slopes and therefore to be 
more centrally concentrated than the satellites of late-type 
hosts. In contrast to this. Sales _& Lambas (2005) using data 
from 2dFGRS, found that the distribution of the satellites 
of red hosts has a shallower profile than the satellites of 
blue hosts, which even deviates from a power l aw, fia t tening 
at small projected separation. More recently IChenI (|2008l ) 
using SDSS data, found that the satellites of both red and 
blue hosts follow similar power slopes after correcting for in- 
terlopers (galaxies mistaken as satellites through projection 
but not actually physically bound to their host). 

Divi ding th e ir sat ellite sample into red and blue pop- 
ulations, ICherJ (|2008l ) also finds that red satellites are 
more centrally concentrated than blue satellites, a trend 
which is also seen in the semi-analytic galaxy samples pro- 
duced bv lSales et al.l ([2007) using the Millennium Simula- 
tion. One explanation of this is that red satellites were ac- 
creted into their host's halo at earlier times than the blue 
satellites, which is consistent with the observational find- 
ing that red satellites have an anisotropic angular distri- 



bution with a preferenc e of b eing aligne d along the ma- 
jor axes of their hosts IIBrain erd 2005; Yang et alj 20061: 
lAzz aro et al.ll2007l ; lBailin et aLll200 8; Agustsson & Brainerd] 
1 20101 ). Comparisons between the radial distribution of satel- 
lites and the dark matter distribution produced from sim - 
ulat ions have also been condu cted by IChen et al.l (|2006r ) 
and Ivan den Bosch et al.l (|2005l ). who find that satellites are 
more centrally concentrated than dark matter sub-halos, but 
consistent with the dark matter profile. 

Satellite galaxies not only provide useful information 
about the formation of galaxies but also about the processes 
that govern galaxy evolution in localised environments on 
scales of ~ 1 Mpc. In the current theory of galaxy evolu- 
tion it is thought that virtually all galaxies start off as blue, 
late-type discs which are then transformed by various pro- 
cesses into red, early types. This is supported by the ob- 
served bimodalit y of galaxies which ca n be seen out to z ~ 1 
(|Bell et al.ll2004l ; IWillmer et al.ll2006l; iPrescott et al.ll2009l) 
and s tudies such as IWillmer et al.l (|2006l 'l and iFaber et al.l 
(|2007|) which have shown there has been a doubling in the 
stellar mass density of galaxies on the red sequence over the 
last 7-8 Gyr. 

The main processes believed to be responsible for trans- 
forming blue/late-type galaxies into red/early types involve 
the quenching of star formation, and include maj or-mergers 
JToomre fc Toomrd 1 19721 ; iHopkins et" aL '2008a^ ), feedback 
from active galactic nuclei ( Bower et al. 2006 ; Croton et al.l 
120061 ) and the depletion of gas reservoirs that fuel star for- 
mation. For satellite galaxies, the dominant process is most 
likely to be gas depletion, caused by the stripping of gas via 
a number of different hydrodynamical and radiative interac- 
tions with their hosts, acting over different timescales. 

When a satellite halo is accreted by the larger halo 
of its host, hot gas from the satellit e may be removed 
in the process kno wn as strangulation (jLarson et al.lll980l ; 
iBalogh et al.ll200d ). resulting in the gradual decline in star 
formation over long timescales (> 1 Gyr), as its fuel for fu- 
ture star formation is depleted. Star formation can be shut 
off more rapidly if the satellite is subjected to sufficient ex- 
ternal pressure that its cold gas r eservoir is removed in the 
process of ram-pressure str ipping (IGunn fc Gottjl 19721) . Gas 
stripping via harassment (jMoore et al.l 119961 ) whereby the 
dark matter sub-halos of satellites are heated after undergo- 
ing frequent high-velocity encounters with other dark mat- 
ter halos, is also a possibility, although this process is more 
likely to occur in galaxy clusters rather than small groups. 

Recent studies using the SDSS have indicated that 
strangulation is the main mechanism causing the transition 
of satellites from the blue to t he red sequence. Producing 
a group catalogue from DR2, IWeinmann et al.l (|2006l ) in- 
vestigate how the fractions of early- and late-type satellites 
vary as a function of halocentric radius, halo mass and lumi- 
nosity, observing that the early-type fraction increases with 
decreasing halocentric radius, increasing halo mass and in- 
creasing luminosity. They argue that the increase in early- 
type fraction with luminosity at fixed halo mass is not ex- 
pected if ram-pressure str ipping or harassm ent is the pri- 
mary cause of gas removal. lAnn et al.l (|2008J ') use SDSS DR5 
to investigate how the early-type fractions of satellites sur- 
rounding isolated hosts vary as a function of luminosity and 
projected distance. They find that the early- type satellite 
fraction increases significantly with decreasing projected ra- 
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dius for early-type hosts and stays approximately constant 
for late types. They conclude that hot X-ray emitting gas of 
the early-type hosts is responsible for the removal of gas. 

Using an S PSS DR4 group catalogue, 
Ivan den Bosch et al.l (|2008l ) compare the concentra- 
tions and colours of centrals and satellites of the same mass. 
By matching central and satellite pairs in both stellar mass 
and concentration, they find there is a significant difference 
in colour. Centrals and satellites matched in both stellar 
mass and colour, on the other hand, show no difference 
in concentration. Under the assumption that centrals are 
the progenitors of satellites (centrals change into satellites 
after being accreted into a larger halo), this implies that 
either the strangulation or ram pressure stripping process is 
occurring rather than harassment, which is believed to have 
a significant effect on the morphology of galaxies. Investi- 
gating the red fraction of satellites by mass they estimate 
that 70 per cent of satellite galaxies with A1* = IO^AIq 
have undergone satellite quenching in order to be on the 
red sequence at the present, with 30 per cent already red 
before becoming a satellite. For more massive satellites 
they find quenching to be less effective, with 65 per cent of 
satellites with M, — IO'^'^A^q being red before accretion 
and virtually all satellites with 7V(, — W^^A4q already 
being red before becoming satellites. 

Furthermore, satellites are believed to affect the evolu- 
tion of their hosts. Minor mergers between dwarf satellites 
and their central hosts p rovide one way in which to distort 
and thicken galaxy discs (JQuinn et al. 199 3) and enlarge the 
bulge components of discs ( Dommgucz-P almero fc Balcellsl 
120081 ) , and mergers involving gas-rich satellites may also pro- 
vide gas which co uld prolong or repleni s h star formation in 
early- type spirals (| White fc Frenklll99ll : [Hau et al.ll2008l ). 

Finally, satellite systems analogous to the Milky Way 
(MW)-Magellanic clouds systems have also become of much 
interest lately due to their apparen t rarity. Observat ional 
studies such as James fc Ivorvl (|201l|) , iLiu et al.l (|201ll ') and 
iTollerud et"ai] l|201ir ) find that onlv ~10 per cent of MW Uke 
hosts have one Magellanic cloud like satellite and only ~5 
per cent have two. Similarly N-body simulations have shown 
that less than 10 per cent of MW sized dark matter haloes 
contain two Magellanic cloud siz ed sub-halos (jBusha et al.l 
l20ld : lBovlan-Kolchin et al.ll2010h . 

In this paper we use data taken from the Galaxy And 
Mass Assembly (GAMA) survey to investigate: firstly the 
radial distribution of satellites surrounding a sample of iso- 
lated host galaxies, as a function of host mass and colour; 
and secondly how the red fraction of satellites depends on 
the projected distance between satellite and host, the stellar 
mass of the satellite and the stellar mass of the host. The 
GAMA redshift survey extends up to 2 mag deeper than the 
SDSS Main Galaxy Sample and, unlike previous studies, we 
use isolation criteria that take account of stellar mass esti- 
mates. The results are discussed in terms of the quenching 
mechanisms that act on satellites. 

The structure of this paper is as follows. In Section 2 we 
outline the GAMA survey from which our samples of galax- 
ies are taken. In Section 3 we define our criteria used to se- 
lect isolated hosts and satellite galaxies. We also define how 
we select a control sample of neighbours. Section 4 describes 
how we divide the red and blue populations of galaxies, com- 
pare properties of the hosts and satellites, and determine the 



projected density of satellites as a function of host mass and 
host colour. We show how the red fraction of satellites de- 
pends on satellite mass, host mass and projected radius and 
discuss potential processes that could produce these results 
in Section 5. Finally in Section 6 we summarise our main re- 
sults. Throughout this paper we assume values of Ho ~ 70 
kms"^ Mpc~\ Qm = 0.3 and Qa = 0.7. 



2 DATA 

2.1 Galaxy And Matter Assembly 

GAMA is a project to construct a multi-wavelength (far-UV 
to radio) database of ^ 375 000 galaxies, by combining pho- 
tometry and spectroscopy from th e latest wide-field survey 
facilities l|Driver et al.ll2009l . IJOllI '). Currently covering 144 
deg^ and going out to z ~ 0.5, GAMA will allow the study of 
galaxies and cosmology on scales between 1 kpc and 1 Mpc 
and provide the link between wi de-shallow surveys, s uch as 
the SDSS Main Galaxy Sample (Stra uss et ai]|2 002''). 2dF- 
GRS (CoUess et al. 2001 ) and 6dF GRS (I Jones et al. 20041 ). 
and na rrow-deep surveys D EEP2 (jDavis et al.l 120031 ) and 
VVDS (|Le Fevre et al.ll2005l ). 

Central to GAMA is a redshift survey conducted at 
the 3.9-m Anglo-Austral ian Telescope (A AT) using the 
AAOmega spectrograph (jSharp et al.l 120061 ). which is cru- 
cial to addressing the main objectives of the project. These 
include determining the dark matter halo mass function of 
groups and cluster s (,Ek c ct al. 2006), measuring the stel- 
lar mass function IjBaldrv et al.l 120081 ) of galaxies down to 
Magellanic Cloud masses and determ ining the recent galaxy 
merger rate (jDe Propris et al.ll2005l ). 

The redshifts used in this paper were obtained as part 
of the initial spectroscopic survey known as GAMA I, car- 
ried out over 66 nights between March 2008 and May 2010. 
This consists of three 12 x 4 deg fields at 9h, 12h and 14. 5h 
(G09, G12 and G15) along the celestial equator and cover- 
ing in total 144 deg^. For detailed descriptions of the spec- 
troscopic target selection and th e tiling strate g y use d for 
GAMA, the reader is referred to iBaldrv et all |2010l ) and 
iRobotham et al.l (|2010l ). respectively. In brief, galaxies are 
selected for spectroscopy using an input catalogue drawn 
from the Sloan Digital Sky Survey (SDSS) Data Release 6 
lAdelman-McCarthv et al. 2008) and UKIRT Infrared Deep 
Skv Survev (UKIDSS) (|Lawrence et al.ll2007l ). 

In the following analysis we use data for galaxies which 
make up the r-band limited Main Survey, which contains in 
total 114 441 galaxies, which are spectroscopically selected 
to have Galactic extinction corrected Petrosian magnitudes 
(|Petrosianlll"97^ of rpetro < 19.4 in fields G09 and G15 and 
rpetro < 19.8 in G12. 

The high density of spectra per square degree of sky and 
high completeness of the redshift surv ey (98 per cent down 
to rpetro = 19.8, [Driver et al.l l201ll ) required to achieve 
the main objectives of the project, make GAMA an ideal 
dataset to study satellite galaxies, since it does not suffer 
from the same incompleteness as other spectroscopic sur- 
veys due to fibre collisions. The SDSS for example has a 
minimum fiber spacing of 55 arcsec, resulting in 10 per cent 
of SDSS targets being missed from the spectroscopic sam- 
ple because each area is generally tiled only once, or twice 
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in the overlap between plates. The fraction missed is higher 
in higher density regions, e.g., for galaxies with two other 
targets within 55 arcsec, the chance of obtaining an SDSS 
redshift is half that of galax ies with no close neighbours 
(figure 3 of lBaldrv et~alll2006l '). This is not true for GAM A 
because each area is tiled 4 or more times, close targets are 
given a higher pri ority in early visits, and high completeness 
is a primary goal (jRobotham et al.ll2010r ). 



2.2 Distances 

From the 114441 galaxies in the r-band magnitude limited 
Main Survey we choose galaxies with reliable redshifts (red- 
shift quality values of Q ^ 3), in the range 0.01 < z < 0.15, 
which results in a sample of 34 102 galaxies, from which we 
search for isolated galaxies and satellites. In the low redshift 
regime of this galaxy sample, the recessional velocities of 
galaxies are significantly affected by peculiar motions which 
can cause distance estimat es to be in error, if s imply assum- 
ing Hubble flow velocities ({Masters et al.ll2004l '). To mitigate 
the effects of peculiar motions on distance we make use of 
a parametric model of t he local v elocity field known as the 
multi-attractor model of iTonrv et al., (,2000. ). 



2.3 Photometry 

The photometry used in this paper includes gri Petrosian 
magnitudes taken from SDSS DR6 and Kron-like AUTO 
magn itudes measured using SExtractor (jBertin fc Arnoutd 
Il99a ). fro m our own re-re duction of the SDSS images de- 
scribed in lHill et al.l (|201in . In brief the photometry for the 
SDSS is obtained for five broad-band filters [ugriz) using 
a dedicated 2.5-m telescope at Apache Point, New Mexico, 



equipped with a mosaic CCD camera ( Gunn et al.lll998l ) 
and calibrated with a 0.5-m telescope l|Hogg et al.ll200lh . 
Fo r greater detail re gardi ng the SDSS the reader is referred 
to I York et al] l|2000l ) and lStoughton et all (|2002l '). 

Unless otherwise stated, absolute magnitudes are cal- 
culated from GAM A AUTO mags such that: 



M = ruAUTO - 5 login Dl-25- A- K 



(1) 



where mAUTO is the apparent AUTO magnitude of a galaxy, 
Dl is the lu minosity distance in M pc, A is the Galactic ex- 
tinction from ISchlegel et al.l (|l998 ) and K the K-correction 
to z = 0.0, determined usi ng the K-Correct v_4_l_4 code 
of iBlanton fc Roweid (|2007l '). For the ~ 700 galaxies where 
SExtractor has failed to assign an AUTO mag in any of 
the g,r,i-hanAs we use the Petrosian magnitude (mpetro) to 
calculate absolute magnitudes. 

In the next sections we describe the criteria and meth- 
ods used in the search for isolated and satellite galaxies. 



3 SELECTION OF SATELLITE SYSTEMS 

3.1 Isolated Galaxies 

Before searching for satellites we must define a set of cri- 
teria to search for isolated galaxies which may host them. 
In previous papers isolated galaxies are usually found by 
searching for galaxies which have no neighbours brighter 



than a given magnitude (either apparent or absolute) con- 
tained within a cylinder, determine d by a velocity and pro- 
jected rad ius (^Zaritsky et al. 1993; ISales fc Lambad [200J; 
iBailin et al. 2008 ; Ann ct al. 2008i). Instead of using an iso- 
lation criteria that depends solely on luminosity, here we 
use a criterion that also takes into account an estimate of 
the galaxies' stellar masses, which makes the selection more 
physically motivated. 

To determine stellar masses we use an expression based 
on the relationship between a galaxy's {g — i) colour and i 
band stellar mass-to- light ratio (Mt/L) from iTavlor et al.l 
(|201ll ). which is given by: 



logioTW. = -0.68 + 0.73(p - i) - OA{Mi - 4.58) 



(2) 



Here A^, is the stellar mass of the galaxy in solar units 
(A1,/A^ solar), {g — i) is the rest-frame colour, and Mi the 
absolute i-band magnitude of the galaxy calculated using 
GAM A AUTO magnitudes. This way of estimating stel- 
lar masses has the advantage of only using two luminosi- 
ties in a transparent way, unlike estimating stellar masses 
determined via SED fitting; and several authors have re- 
cently su^gested__th^A;(*/Lcorrela±es reliably with 
g - i (|Gallazzi fc Bellll2009l : IZibetti et aT]|2009l : iTavlor et all 
l2010t ). 

The central 95 per cent of the rest-frame {g — i) colour 
distribution of the galaxy sample considered in this study 
corresponds to 0.26 < {g — i) < 1.24, implying that the 
galaxy stellar mass-to-light ratio in the i-band typically 
varies by as much as a factor of 5 from 0.32 to 1.68. In order 
to be robust against colour errors we restrict the mass-to- 
light ratio to this range. 

In the search for isolated centrals we limit ourselves to 
using galaxies with rp^tro ^ 19.4 in GAMA fields G09 and 
G15, and rpetro ^ 19.8 in GAMA field G12. In order to be 
considered isolated, the centrals must not have any compa- 
rably massive neighbours within a large surrounding region. 
We define isolated galaxies as those which which have: 

(i) A stellar mass more than 3 times that of any neigh- 
bours {Mtjao > SMt, Neighbour), wlthiu a projected radius 
of Rp ^ IMpc and |cAz| sj 500 kms~\ 

(ii) An apparent r-band Petrosian magnitude such that: 



rpetro < rii. 



2.5(logio 3) - Ar 



(3) 



with riim being 19.4 for G09 and G12 or 19.8 for G15 and 
where Ar is given by: 



Ar = 2.5[logi(,(X*/Lr)™^ - \og-^^^{M*/Lr)] 



(4) 



and {M,/Lr)max is a 'maximum' mass-to-light ratio in the 
r-band. r-band mass to light ratios are determined by di- 
viding the stellar masses by absolute r-band Petrosian lumi- 
nosities in solar units. This maximum stellar mass ratio is 
chosen to have the value {Mf/Lr)max ~ 2.14. The purpose 
of this condition is that it ensures that isolated galaxies are 
sufficiently brighter than the limit of the field, such that even 
neighbouring galaxies with the highest mass-to-light ratios 
that are fainter than the limit still have a stellar mass that 
is a factor 3 less than that of the central. 

(iii) A projected distance greater than 0.5 Mpc from each 
edge of the survey regions. This third condition ensures > 
80% of the area of the 1 Mpc circle surrounding a galaxy is 
within the survey region. 
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We find that 7 288 galaxies out of the 34 102 galax- 
ies satisfy the second and third conditions, and overall we 
find 3 536 galaxies to be isolated. As this study is focused 
on satellites hosted by typical galaxies and not those in 
large groups and clusters, we remove all 22 galaxies with 
logj^g A1» ^11.5 that may be brightest cluster galaxies, re- 
sulting in a sample of 3 514 isolated galaxies. 



3.2 Satellite Galaxies 

After finding isolated systems we search for satellites around 
these galaxies. We define satellite galaxies as those sur- 
rounding isolated galaxies, which have a stellar mass that 
is at most one-third that of the central galaxy {A4, < 
1/3A^*, cen), within a projected radius of Rproj ^ 500 kpc 
and \cAz\ ^ 500 km s~^ . We chose a minimum h ost to satel- 
lite mass ratio of 3:1 as I Hopkins et al.l l|2008bl ) find that a 
3:1 mass ratio merger event is the limit for an L* disk to 
survive as a disk gala:xy. A satellite with one-third the mass 
of its host will be ~ 1.2 mags fainter, assuming it has the 
same stellar mass-to-light ratio. 

Out of the 3 514 isolated galaxies we find that 1426 
host a total of 2 998 satellite galaxies. Noting that the more 
massive isolated galaxies are more likely to have satellites, 
and the more massive hosts are more likely to have multiple 
satellites, simply because of the selection, we find that most 
(59.4 per cent) have no satellites and the mean number of 
satellites per isolated galaxy is 0.85. Excluding the isolated 
galaxies which do not host satellites we find that the mean 
number of satellites per central host is 2.10. Varying the 
projected radii and velocity differences in the isolation and 
satellite search criteria results in slightly different numbers 
of satellites per host and per isolated galaxy. 



3.3 Other Small Neighbours 

For the purposes of determining how the properties of satel- 
lite galaxies depend on the host properties, we produce a 
comparative sample of smaller neighbouring galaxies with 
similar masses to the satellites, which satisfy the same crite- 
ria but with projected distances in the range 0.5 ^ Rproj ^ 
1 Mpc. To ensure these neighbouring galaxies have no host, 
we check to see that these neighbouring galaxies have no 
nearby galaxies which are greater than 3 times the stellar 
mass of the small neighbour, within Rp < 500 kpc and with 
jcA^I < 500 kms~^. As these galaxies are selected in a sim- 
ilar way to the satellites, we consider this to be a control 
sample. Surrounding the 3 514 isolated galaxies we find a 
total of 2 304 of these small neighbours. 



4 RED AND BLUE POPULATIONS OF 
GALAXIES 

In this paper we ma k e use of th e well know n colour bimodal- 
ity of galaxies f|Strateva et akl lJ^OOl: Bal drv et al.l [2004') to 
divide the hosts and satellite galaxies into red and blue pop- 
ulations, with a colour-mass diagram (CMD). This enables 
us to compare the properties between the populations, to 
examine how the radial distribution of satellites depends on 




10 
Log(Stellar Mass) 



Figure 1. g-r colour-mass diagram for the entire sample of 34 102 
galaxies within 0.01 < z < 0.15 from the r-band limited Main Sur- 
vey. The data points are weighted by 1/Vmax and represented as a 
logarithmic contour plot. Red and blue populations are separated 
by the line (g - r) = 0.03(logiQ M*) + 0.35. 



host colour, and to determine how the red fraction of satel- 
lite galaxies varies as a function of both projected radial 
distance from the host and stellar mass. 



4.1 Colour- Mass Distributions 

Using g — r calculated from the AUTO mags, we produce a 
CMD using the sample of 34 102 galaxies in the range 0.01 < 
z < 0.15. Figure[l]shows the CMD plotted as logarithmically 
spaced contours. Each data point is weighted by 1/Vmax, 
where Vmax is the maximum comoving volume, within which 
the galaxy co uld lie dependin g on its redshift and the limits 
of the survey (jSchmidtlligli i. 

Colour bimodality can clearly be seen in Figure [TJ and 
we choose to separate the red and blue populations using a 
straight line with the equation: 



(g - r) = 0.03(logio M,) + 0.35. 



(5) 



Figure [T] lo oks almost identical to the CMD in figure Al of 
Ivan den Bo sch et al. (2008), produced using stellar masses 
derived from lBell et all (|2003i ). CMDs for the isolated galax- 
ies, central hosts and satellites with projected distances of 
Rproj ^ 500 kpc can be seen in Figure[21 which clearly shows 
that the majority of isolated galaxies which host satellites 
are red (906 out of 1 426 hosts). Out of the 2 998 satellites we 
find 1 097 red and 1 901 blue galaxies. By further dividing 
the sample into red/blue satellite/host systems, we find that 

41.2 per cent (1 235) of the satellites are blue with red hosts, 

33.3 per cent (998) are red satellites with red hosts, 22.2 per 
cent (666) are blue satellites with blue hosts, and only 3.3 
per cent (99) are red satellites with blue hosts. Overall the 
blue fraction of satellites is 63.4 per cent, which rises to 87.1 
per cent for blue hosts compared to 55.3 per cent for red 
hosts. 
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Figure 2. g-r colour-mass diagrams for (a) the 3 514 isolated 
galaxies (b) the 1 426 central hosts and (c) 2 998 satellites. 
Red and blue populations are separated by the line {g — r) = 
0.03(logio A^,) + 0.35. 



4.2 Luminosity and Mass Distributions 

In this section we compare the properties of the different 
samples of galaxies. In Figure[3]we show the distributions in 
himinosity and stellar mass for red and blue hosts and satel- 
lites. This shows, as expected, the red population on average 
being more luminous and massive than the blue popula- 
tion. We calculate mean stellar masses and absolute r-band 
magnitudes for red host galax;ies as logjQ M* — 10.81 and 
Air = —21.81 respectively, compared to logjQ A4* = 10.31 
and Mr — —21.28 for the blue hosts. The latter is sim- 
ilar to the value for the Milky Way, and the Schechter 
break (M*) of blue galaxies (Loveday et al. 2011 in prepara- 
tion). Similarly for the red satellites we calculate means of 
logioX. = 10.00 and Mr = -19.96, and log^, X. = 9.35 
and Mr = —19.30 for the blue. Both populations are typi- 
cally more massive and luminous than the Large Magellanic 
Cloud. 

Comparing the host galaxies to isolated galaxies with- 
out satellites, we find that both the blue and red galaxies 
with satellites are more massive and more luminous than 
those without, which is most likely a result of our satellite 
selection as the more massive isolated galaxies have a larger 
chance of hosting a satellite. The average mass and luminos- 
ity of the red galaxies without satellites is logj^Q A4, — 10.61 
and Mr = —21.39 compared to the blue galaxies which have 
logioX, = 10.12 and Mr = -20.94. 

The mean masses and luminosities of the small neigh- 
bours are very similar to those of the satellites and we find 
means of logjg A^, = 10.77 and Mr = —20.17 for red galax- 
ies and logioA^. = 9.42 and M,. = -19.48 for the blue. 
This is reassuring as it means that the small neighbours do 
indeed have similar properties to the satellites and can be 
considered a 'control' sample, allowing comparisons between 
the two. 

Other quantities useful for comparison with cosmolog- 
ical models include the stellar mass and luminosity differ- 
ences between hosts and satellites. In Figure |4] we show his- 
tograms of these quantities for the blue and red hosts. 

Due to our satellites having stellar masses that are less 
than 1/3 of their host mass, AlogjpA^, has a minimum 
value of 0.477, and we find a maximum of Alogj^pA^* = 
3 (1/lOOOth the mass of the host). We find find that the 
average satellite in this study is ~ 1/lOth of the mass of the 
host and calculate mean logarithmic mass differences of 1.05 
and 1.39 for red and blue satellites with red hosts, and 0.95 
and 1.11 for red and blue satellites with blue hosts. 

As for luminosity difi'erence, we find that satellites are 
typically ~ 2 magnitudes fainter than their hosts. Red and 
blue satellites of red hosts, are both found to have mean 
magnitude differences of 2.37, and red and blue satellites 
with blue hosts have magnitude differences of 2.54 and 2.21. 
We also find a small fraction of satellites which have a neg- 
ative value of AMr and are brighter than their hosts, which 
would be missed by other studies. In other words, when con- 
sidering red hosts, they can be more than three times more 
massive than their blue satellites yet be of similar luminos- 
ity. 
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Figure 3. Histograms showing the distribution of r-band absolute magnitudes (a,b) and stellar masses (c,d) for the blue and red 
populations of satellites with Rproj ^ 500 kpc (solid lines) and host galaxies (dashed lines). 
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Figure 5. (a): Projected densities of satellites with Rproj ^ 
500 kpc, divided into four host stellar mass bins, using projected 
radius bins of width AR = 0.05 Mpc. (b); The same as above 
corrected for interlopers. The error bars shown are Poissonian. 



5 PROJECTED DENSITY OF SATELLITES 

To investigate the radial distribution of the satellites we cal- 
culate the projected density per host given by: 

Nsat{R) 



EiR) = 



1^NHost{R2^ - Rl) 



(6) 



where Nsat{R) is the number of satellites found within a 
shell bounded by inner and outer projected radii, Ri and 
7?2, from the host and Nnost is the number of host galaxies 
and R is the midpoint of the bin. 

Using projected radius bins of width AR = 0.05 Mpc 
and dividing the satellites into four host stellar mass bins; 
9.5 < logioX* < 10.0, 10.0 < logioX, < 10.5, 10.5 < 
log^^M, < 11.0 and 11.0 < log^QM* < 11.5, we produce 
Figure[5ja) which shows how the projected density of galax- 
ies varies as a function of projected distance. We fit the pro- 
jected density with a power law, parametrized by a slope, 
a, and normalisation A: 

E(i?) = yl(i?/Mpc)°. (7) 

with A in Mpc~^ as S is in Mpc"'^. 

In Table [T] we show the best-fitting parameters and 
errors for each of the curves, and find that there is little 
variation in the best-fitting slope, which steepens slightly 
from a = —0.94 ± 0.04 for satellites with hosts in the 




Figure 6. Projected densities of satellites with Rproj ^ 500 
kpc, for the red and blue hosts with masses in the range 9.5 ^ 
logjg Mt ^ 11.5 , using projected radius bins of width AR = 0.05 
Mpc. The interloper corrected values are connected by dotted 
lines. The error bars are Poissonian. 



mass range 11.0 < log^pAl, < 11.5, through to a = 
— 1.06 ± 0.03 for satellites with hosts in the mass range 
9.5 < logjgA^, < 10.0. Dividing the satellites into those 



with red and blue hosts (Figure |6} , with host masses be- 
tween 9.5 < logj^Q A4f < 11.5, we find a best fitting slope of 
a — —0.94 ± 0.03 for satellites with red hosts compared to 
those with blue host which have a = —1.04 ± 0.04. These 
slopes suggest that the blue, low mass host galaxies have 
slightly more centrally concentrated satellites. The slopes 
however are affected by interlopers, galaxies which are mis- 
taken as satellites by projection and not actually bound to 
their hosts. 

In order to correct for interlopers we use a simple 
method of subtracting the average projected density of the 
small neighbours, with projected radii in the range 0.5 < 
Rproj < 1.0 Mpc, from the projected density of the satel- 
lites in each bin: 



E(i?)i 



E(i?)sat - (E(0.5 <R< 1.0)emallncigh> • (8) 



Note that this assumes that these galaxies in the 0.5-1 Mpc 
annuli represent the average large-scale or 'super cluster' 
environment around the hosts, and the correction will be 
an overestimate in hosts where a significant bound satellite 
population extends beyond 0.5 Mpc, in particular for our 
most massive host sample (11.0-11.5). 

The result of applying this interloper correction and fit- 
ting the slope for each of the host mass ranges (and colour) 
can also be seen in Table[T] which shows that the best-fitting 
slope now significantly varies with host mass and steepens 
from a — —1.09 ± 0.04 for satellites with hosts in the mass 
range 11.0 < logj^pAl* < 11.5 through to a = — 1.4±0.1 for 
satellites with hosts in the mass range 9.5 < logj^Q M < 10.0. 
For the red and blue hosts, we obtain interloper corrected 
slopes of a = —1.13 ± 0.03 and a — —1.38 ± 0.04 respec- 
tively. The results of this inteloper subtraction can be seen 
in Figure [5jb) and as the dotted lines in Figure [G] 

An observed decline in the projected density, which 
overall can be approximated by S(7?) oc R~^ , is simi- 
lar to other studies that inv estigate the radial di stribu- 
tion of satellites. For example, iLorrimer et al.l ([199J) using 



© 2010 RAS, MNRAS OOO.nHTSl 



GAM A: Satellite Galaxies 9 

Table 1. Best fitting A and a parameters (Equation [7| for the projected densities in Figure [5] and Figure |6l corrected and uncorrected 
for interlopers. The uncertainties quoted are derived from least-squares straight-line fits to log S versus log R with the data points shown 
in the figures. 



Host Stellar Mass Range A^Sattot Auncorr (Mpc ^) 



At, 



(Mpc-2) 



aintc 



9.5 ^ 



,A1. s£ 10.0 



10.0 ^ logio M, ^ 10.5 
10.5 ^ logio A4, ^ 11.0 
11.0 ^ logioA4* ^ 11.5 



111 


0.14 ±0.03 


457 


0.13 ±0.01 


1225 


0.20 ±0.01 


1140 


0.23 ±0.02 


700 


0.14 ±0.01 


2 223 


0.21 ±0.01 



-1.06 ±0.10 
-1.09 ±0.05 
-0.94 ±0.03 
-0.94 ±0.04 



0.05 ±0.01 
0.04 ± 0.04 
0.10 ±0.01 
0.14 ±0.01 



-1.40 ±0.10 

-1.46 ±0.05 
-1.14 ±0.04 
-1.09 ±0.04 



ALL BLUE HOSTS 
ALL RED HOSTS 



-1.04 ±0.04 
-0.94 ±0.03 



0.05 ±0.01 
0.12 ±0.01 



-1.38 ±0.04 
-1.13 ±0.03 



the Centre for Astrophysics redshift survey, find a slope of 
a — —0.91 ±0.05 for satellites with magnitudes in the range 

— 16 < B < —18, after correcting for background galax- 
ies, which is slightly shallower than the interloper corrected 
slopes we obtain. Furthermore, they divide their satellite 
sample into early- and late-type hosts, to find that satellites 
with early-type hosts are more centrally concentrated than 
late types, finding slopes of a = —1.01 ± 0.1 for satellites of 
early-type hosts and a — —0.78 ± 0.05 and for the late-type 
hosts, which is the opposite to what we find. Although our 
slope for satellites of red hosts is consistent with their re- 
sult for early-type hosts, our result for the satellites of blue 
galaxies differs from their findings for late-type hosts. 

More recentlv. lChen et al.l (|2006l ) find an interloper cor- 
rected slope of a = —1.58 ± 0.11, using a galaxy sample 
from SDSS DR4. Here satellites are defined as galaxies with 
Air < —18.0, which are 2 mags fainter than their hosts, 
with Rproj =S! 500 kpc and AV ^ 500 kms~^ and isolated 
galaxies are defined as those with M,- ^ —20.0, have no 
nearby galaxies within Am ^ 2 mag, AV ^ 1000 kms^^ 
and Rproj ^ 500 kpc. Using these same sele ction c riteria 
and a larger data sample from SDSS DR6, IChenI (|2008l ') 
find that the interloper corrected slopes of satellites of red 
and blue hosts are almost identical and obtain values of 
a = -1.54 ± 0.12 and a = -1.50 ± 0.12 for satellites of blue 
hosts and red hosts, respectively. Although these values are 
similar to our slopes for the blue and low-mass hosts , they 
are qu ite different for the red and high-mass hosts. IChenI 
l|2008r) and others have used r-band luminosity for isolation 
criteria. For red hosts, this is significantly more stringent 
than using an estimated stellar mass that allows for equally 
luminous blue satellites. This could explain some of these 
differences. 

I Sales et al.l l|2007l ) using a sample of galaxies from the 
semi-analytic Millennium Simulation find a slope of a = 

— 1.55 ± 0.08, for satellites again that are at least 2 magni- 
tudes fainter than their hosts, within Rproj ^ 500 kpc and 
AT/ ±500 kms-1. 

Investigating t he satellites of 34 elliptical galaxies 
iMadore et al.l (|2004l ) find a decline with a slope of a = —0.5 
out to Rproj — 150 kpc, which flattens out at greater pro- 
jected radii. Although the slope they obtain is again incon- 
sistent with our results, this may be due to the small num- 
bers of galaxies involved. They do note that there is a flat- 
tening in the projected density at Rproj > 300 kpc, which 
can also be seen in our data. 



Finally, lAnn et al.l (|2008l ') using SDSS DR5 data, de- 
flne satellites as galaxies with Mr < —18.0, which are more 
than 1 mag fainter in the r-band than their hosts, within 
Rproj < 1 Mpc and AV sj 1000 kms~\ They also use 
an isolation criterion, so that isolated galaxies are required 
to have Mr < —19.0, and no neighbours within Am ^ 1 
mag of the host, within a Rproj ^ 500 kpc and AV ± 500 
kms~^, producing a sample of 2 254 hosts with 4 986 satel- 
lites. Without correcting for interlopers they flt slopes for 
satellites within Rproj ~ 200 kpc of their hosts, and ob- 
tain values of Q ~ —1.8 for both the late- and early-type 
satellites of late-type hosts and for the early-type satellites 
of early-type hosts. For the late-type satellites of early-type 
hosts a shallower slope of a ~ —1.5 is found. 

The shallower slopes for the high-mass hosts com- 
pared to low-mass hosts measured over the range 0.05- 
0.5 Mpc are expected i f satellite density proflles are r elated 
to dark-matter proflles (jNavarro. Frenk fc Whitell996l. here- 
after NFW; surface density proflle given by iBartelmannI 
[l99a). The higher- mass hosts would be expected to have 
larger NFW scale radii and therefore shallower dark-matter 
proflles over the radial range probed by our study. 



6 RED FRACTION OF SATELLITE GALAXIES 

After measuring the projected density of satellites, we de- 
termine how the red fraction of satellites varies as a function 
of mass and projected distance from the hosts. It should be 
noted that the most robust results, free of signiflcant selec- 
tion effects, are the differences between the satellites and 
small neighbours. 



6.1 Red Fraction as a Function of Projected 
Distance 

To measure the red fraction of satellites {Rproj ^ 0.5 Mpc) 
and small neighbours (0.5 ^ Rproj ^ 1 Mpc) we divide the 
number counts of red galaxies by the total number of galax- 
ies in bins of projected radius of width ARproj =0.1 Mpc. 
By dividing the sample of satellites into those with red and 
blue hosts we produce Figure EJa). By dividing the sam- 
ple of satellites into four bins of host stellar mass as before, 
we produce Figure [TJb). The error bars on the red fraction 
of galaxies in all plots are la beta distr i bution confldence 
intervals explained in detail in ICameronI (|2010|). These are 
an improvement over other methods for estimating binomial 
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Figure 7. (a): The red fraction of satellite galaxies as a func- 
tion of projected radius for the red and blue hosts, (b): The red 
fraction of satellite galaxies as a function of projected radius, for 
four different host stellar mass ranges. We define galaxies with 
Rproj ^ 0-5 Mpc as satellites and galaxies with 0.5 ^ Rproj ^ 
1.0 Mpc as small neighbouring galaxies which have been selected 
in a similar way. The error bars shown are la beta distribution 
confidence intervals. 



confidence intervals, such as Poisson errors, which are liable 
to misrepresent the degree of statistical uncertainty for small 
samples, or fractions that approach values of or 1. 

From Figure [TJa) it is clear that the red fraction of 
satellites of red hosts is significantly higher than those with 
blue hosts. The average red fraction of all the points between 
0.05 ^ Rproj ^ 1 Mpc for the satellites of red hosts is 0.38 
compared to 0.13 for the blue hosts. The red fraction of 
satellites for red hosts increases by a factor of 1.49 from 0.37 
at Rpro] = 0.45 Mpc to 0.55 at 0.05 Mpc. For the satellites 
of blue hosts no trend in the red fraction as a function of 
radius is seen. 

Similar trends are refiected in Figure [TJb) as expected. 
The red fraction of satellites and small neighbours increases 
significantly as a function of host mass, from an average 
red fraction of all the points between 0.05 ^ Rproj ^ 



1 Mpc of 0.09 
9.5 s; login X. 
10.0 < logioX. 
10.5 ^ logioAI. 
hosts with 11.0 



for galaxies associated with hosts with 

^ 10.0, 0.08 for galaxies of hosts with 

, ^ 10.5, 0.27 for galaxies of hosts with 

, ^ 11.0, through to 0.54 for galaxies of 

^ logipAl, ^ 11.5. The highest two host 



mass ranges show that an upturn in the red fraction oc- 
curs within Rproj ~ 500 kpc. For the highest host mass 
bin 11.0 ^ logjQ Al ^ 11.5, the red fraction of satellites in- 
creases by a factor of 1.4 from 0.54 at Rproj ~ 0.45 Mpc 
to 0.74 at 0.05 Mpc and for satellites with host masses of 
10.5 ^ logiQ M* ^ 11.0, with the red fraction increasing by 
a factor of 1.7 fi-om 0.22 at Rproj = 0.45 Mpc to 0.37 at 0.05 
Mpc. No trend in the red fraction as a function of radius is 
observed for the 2 lowest host mass bins. 

Following on from this we investigate how the red frac- 
tion varies as a function of satellite mass, by subdividing 
galaxies in the host mass bins above further into satellite 
mass bins. It is important to do this in order to verify that 
the satellite trends are not simply due to mass segrega- 
tion, and to separate 'environ ment quenching' from 'mass 
quenching' (IPeng et al.l |2010|). Figure [8ja,c,e) shows how 
the red fraction varies as a function of radius, with the 
data separated into 3 different plots for three different host 
mass bins, and then subdivided into satellite mass bins of 
9.5 < logipX. < 10.0, 10.0 < logioA^. < 10.5, 10.5 ^ 
logigX* ^ 11.0 and 11.0 ^ logj^^M* < 11.5, depending 
on the host mass bin. The mean numbers of satellites and 
small neighbours in each bin are 36.2, 65.6 and 48.5 for those 
with host masses in the ranges 10.0 ^ log^gA^* ^ 10.5, 
10.5 < logioX. < 11-0 and 11.0 < logmAl* < 11-5 re- 
spectively. In Table [2] we show the average red fractions of 
satellites and small neighbours for each of the different host 
and satellite mass ranges. From these it can be seen that: 

• For satellites with hosts in the range 10.0 ^ log^g Mt ^ 
10.5, little change in the red fraction as a function of radius 
is observed, with the red fraction of the satellites and small 
neighbours remaining around the 0.05 level for galaxies with 
9.0 ^ logigA^* ^ 9.5 and ~ 0.14 for satellites with stellar 
masses in the range 9.5 ^ log^o M, ^ 10.0. 

• For satellites with hosts in the range 10.5 ^ log^g -M^ ^ 
11.0, the red fraction increases with satellite mass, such that 
the average red fraction increases from 0.11 ±0.03 for satel- 
lites with masses 9.0 sj log^jX. ^ 9.5, to 0.33 ± 0.03 
for satellites with masses of 9.5 ^ log^oA^* ^ 10.0, and 
0.57 ± 0.02 for satellites with masses of 10.0 ^ logjo -Mt ^ 
10.5. It also becomes the apparent that the red fraction in- 
creases as distance decreases. The average red fraction of 
the satellites is greater than the average red fraction of the 
small neighbours in the same mass range. For galaxies with 
9.0 < logio M* s; 9.5 this diflterence is 0.05±0.04, 0.15±0.04 
for galaxies with 9.5 < log^Q M, < 10.0 and 0.06 ± 0.04 for 
galaxies with 10.0 < log^gM* ^ 10.5. 

• Satellites with host masses 11.0 ^ logmA^* $; 11.5 
show similar trends and the radial dependence becomes even 
more apparent, especially for the low mass satellites. The 
average red fraction of satellites increases from 0.37 ± 0.07 
for satellites with 9.0 sj log^^M, < 9.5, through to 0.90 ± 
0.02 for the satellites with 10.5 ^ logioX* ^ 11.0. The 
difference between the red fraction of satellites and the small 
neighbours is observed to be 0.25 ± 0.07 for galaxies in the 
mass range 9.0 ^ log^Q A1« ^ 9.5, 0.18 ± 0.07 for galaxies 
with 9.5 < login TW. ^ 10-0. 0-15 ± 0.05 for galaxies with 
10.0 < logic X. =S 10-5 and 0.15 ± 0.04 for galaxies with 
10.5 ^logipX, < 11.0. 

Our results here support the 'galactic conformity' phe- 
nomenon, whereby red/early-type centrals have a signifi- 
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Figure 8. Left panels (a,c,e) show the red fraction of satellite galaxies as a function of projected radius and satellite stellar mass, shown 
for 3 different central host stellar mass ranges. Host masses of 10.0 C logio -M* ^ 10.5, 10.5 ^ logio J^* ^ H-O ^^id H-O ^ l°gio -^* ^ 
11.5 with satellites of stellar masses 9.0 ^ logic -^* ^ ^-^ (''^'^ stars), 9.5 ^ logio -^* ^ l^^-O (green diamonds), 10.0 ^ logio -^* ^ 1^0.5 
(blue triangles) and 10.5 S^ logio -^* ^ ^^-0 (purple squares). Right panels (b,d,f) show the red fraction of satellites with Rp^oj ^ 0.5 
Mpc and small neighbours with 0.5 ^ Rproj ^ 1-0 Mpc as a function of satellite mass, for hosts in the same mass ranges opposite. The 
errors shown are la beta distribution confidence intervals. 



ca ntly hisher fraction of re d/early-type satellites, noted first 
by IWeinmann et aU (|2006|) and also observed by lAnn et aU 
l|2008r ). IWeinniann et al.l (|2006l ) produce a large catalogue of 
~ 53 000 groups, from a sample of ~ 92 000 galax ies from the 
SPSS DR2, using the halo-based group finder of lYang et al.l 
l|2005l) . Defining centrals as the brightest group members 
and satellites as the remaining group members, they deter- 



mine the projected halocentric distances (normalised to the 
virial radius of the system Rvir) of the satellites and halo 
masses of the systems. They find that the early-type frac- 
tions of satellites increase as a function of increasing lumi- 
nosity, increasing group halo mass and decreasing halocen- 
tric distance. For halo masses with 14 < logio -^ < ^^i the 
early-type fraction of satellites increases from ~ 40 per cent 
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Table 2. Average red fraction of satellites and small neighbours of galaxies in Figure [S] The error quoted is the standard error. 



Host Mass Range 


Satellite Mass Range 


Mean Satellite Red Fraction 


Mean Small Neighbour Red Fraction 


lO.OsglogioAI. s£ 10.5 


g.OsglogioAI. ^9.5 
9.5 sglogioM* s; 10.0 


0.08 ±0.02 
0.14 ±0.03 


0.04 ±0.01 
0.14 ±0.02 


10.5^1ogioA^* ^ 11-0 


9.0^ logic A^. ^9-5 
9.5^ logic 7K. «S 10-0 
10.0 < logic >1* ^ 10.5 


0.11 ±0.03 
0.33 ±0.03 
0.57 ±0.02 


0.06 ±0.02 
0.18 ±0.04 
0.51 ±0.02 



11.0 s£ logic -M* =£ 11-5 9.0 s£ logic M* ^ 9.5 
9.5 < logic X* «S 10-0 
10.0 < logical* ^ 10.5 
10.5 s£ logic X. ^ 11-0 



0.37 ±0.07 
0.53 ±0.03 
0.79 ±0.03 
0.90 ±0.02 



0.12 ±0.02 
0.35 ±0.05 
0.64 ±0.03 
0.75 ±0.04 



at 0.9R/Rvir to ~ 60 per cent at O.IR/Rvir- For lower-mass 
halos with 12 < logic ■^ < ^3, the early-type fraction of 
satellites is approximately constant at '^ 25 per cent. This 

is a similar picture to what we see in Figure [T] 

In a study more directly comparable to ours. lAnn et al.l 
(|2008l ) investigate satellite systems of isolated hosts, who 
determine the morphologies of their satellite and host sam- 
ples to investigate the early-type fraction of satellite galax- 
ies as a function of radius. Splitting their hosts into early 
and late types, they find that the early-type satellite frac- 
tion for early-type hosts increases from ~ 0.2 at Rproj = 1 
Mpc to ~ 0.6 at Rproj = 0.05 kpc, whereas the early-type 
fraction for late-type hosts remains roughly constant, ~ 0.2 
for all projected radii, which is almost identical to what we 
find in Figure [^ a). By dividing their sample into hosts with 
Mr < -20.5 and -21.0 < Mr < -19.0, they show in their 
figure 2 that the brighter early-type hosts have an early- 
type satellite fraction which is 0.05-0.20 greater in each pro- 
jected radius bin compared to those with fainter hosts. For 
the brighter host bin, the early-type satellite fraction in- 



creases from ~ 20 per cent at Rpr 



1 Mpc to r^ 70 per 



cent at Rproj = 0.05 Mpc, whereas the early-type satellite 
fraction for fainter hosts increases from ~ 20 per cent at 
Rproj ~ 1 Mpc to ~ 60 per cent at Rproj = 0.05 Mpc. 
For the late-type hosts, the early-type satellite fraction is 
roughly constant with projected radius at ~ 25 per cent for 
hosts with A'lr < —20.5 compared to being roughly constant 
at ~ 15 per cent for hosts with —21.0 < Mr < —19.0. Again 
this trend is similar to our result that the red fraction of 
satellites increases with increasing host mass. 



6.2 Red Fraction as a Function of Mass 

An alternative way of showing how the red fraction varies 
as a function of radius can be seen in Figure [S{b,d,f), which 
shows the red fraction of satellites and small neighbours in 
mass bins of width A logic -^ ~ 0.25, for the same host mass 
ranges as on the left-hand side. As expected the red fraction 
of satellites and small neighbours increases with increasing 
stellar mass. Figure [8jb,d,f) clearly shows that host mass 
has the biggest effect on the red fraction of satellites relative 
to the small neighbours. For the satellites with the lowest 
host masses (10.0 ^ logic A^* ^ 10.5), there is virtually no 
difference between the red fraction of the satellites and small 
neighbours, confirming that the red fraction has no radial 
dependence. For the highest mass hosts, the red fraction of 



the satellites is about 0.2 greater than the red fraction of 

small neighbours. 

Usi ng the SPSS DR4 g r oup c atalogue of lYang et al.l 

((20071), Ivan den Bosch et al.l l|2008r ) define the most mas- 
sive members of groups as centrals and the remaining group 
members as satellites, producing a sample of 218 103 centrals 
and 59 982 satellites. They then investigate how the red frac- 
tions of hosts (centrals) and satellites by mass compare. In 
their figure 8, they show that the difference in the red frac- 
tion of centrals and satellites is the greatest at low masses 
(the satellites have red fractions which are 20-40 per cent 
higher for logic A^* ^ 10.0), converging at higher masses. 
The significant red fraction at the lowest satellite masses 
is discrepant with what we observe, although this may be 
explained by the different definitions of satellite galaxies. 



6.3 Discussion 

We have shown that the red fraction of satellites surround- 
ing isolated hosts increases as a function of host stellar mass, 
satellite mass and decreasing projected separation, indicat- 
ing the quenching of star formation in satellite galaxies is 
more efficient in most massive systems and acts within ~ 
500 kpc of the host. This quenching of star formation is 
likely to be due to the removal of gas by a combination of 
radiative and gravitational effects of the host, which have 
been investigated in previous studies. 

One way of shutting off the star formation in 
satellite galaxies is by strangulation IjLarson et al.l Il98d : 
iBalogh et al.ll2000f ). which is thought to occur when a satel- 
lite's halo is accreted into the larger halo of its host removing 
the satellite's hot, diffuse gas and thus causing a slow decline 
(on timescales r > 1 Gyr) in the star formation as its supply 
of c old gas is suppressed, and t urning spirals into SO galax- 
ies. iKawata fc Mulchaevl (|2008l ) simulate the strangulation 
process, showing that it acts on satellites within R < 500 
kpc on timescales of r ~ 2.0 Gyr. They conclude the pro- 
cess is more effective in relaxed groups of galaxies, which 
are dominated by elliptical galaxies surrounded by an X-ray 
emitting intergalactic medium (IGM). Many of the systems 
with high mass hosts found in this study could be described 
as the system above, making strangulation a likely expla- 
nation for our results. This scenario is also consistent with 
the host mass having the largest effect on the red fraction, 
rather than satellite mass. 

Compared to the slow quenching of star formation 
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which occurs in strangulation, rapid quenching could oc- 
cur if enough external pressure is exerted on the satellite 
to re move its cold gas, in the process of ram pressure strip- 
ping (|Gunn fc Gotj Il972l : [Hester 2006 ). Thought to turn 
spiral galaxies into SO types, ram pressure stripping is most 
effective in clusters and groups of galaxies, where the den- 
sity of the hot IGM and velocities of satellites creates suffi- 
cient pressure to overcome the restoring force keeping cold 
gas in galaxies. i Hesteil l)2006l ) finds, using simulations of 
groups, that low mass satellites are more severely stripped 
over larger distances, than larger spiral galaxies which are 
only stripped at small distances from the host. This is not 
what we see in this study with a higher red fraction over 0- 
0.4 Mpc, compared to the small neighbours, in the highest- 
mass hosts for all satellite masses [see Figure [8je)]. 

Tidal stripping of gas from satellites due to the gravi- 
tational effects of the host is likely to occur simultaneously 
with ram- pressure stripping. iMaver et al.l (|2006l ). using N- 
body simulations of Milky Way sized halos surrounded by 
gas rich dwarfs, show that a combination of tidal and ram 
pressure stripping is more effective at removing gas than ei- 
ther process alone. Tidal interactions between satellites and 
their hosts may explain the streams of HI gas observed lo- 
cally as the Magella nic Stream and High Velocity Clouds 
(|Putman et al.ll2003l) and as the extraga lactic HI clouds ob- 
servcd around M31 (JThilker et aLllpOO-J). M33 ( ICrossi et al.1 
[2008 ) and the M81/M82 group dChvnoweth et al. 2001^ 
Shutting off of star formation might be expected to scale 
with the ratio of satellite to host mass, which is not evident 
in our data. Consider satellites that are ~ 0.01-0.1 times the 
mass of the host in Figure [8] [8.5-10 in Fig. ^h), 9-10.5 in 
Fig.Od), 9.5-11 in Fig.IHf)], it is clear that the difference 
between the satellites and small neighbours depends on the 
host mass but not the ratio between the satellite and host 
mass. This is in agreement with the idea that the probabil- 
ity of enviro nment quenching is independent of the satellite 
stellar mass ijPeng et al.ll201ll ). 

Another process which is tho ught to be able to strip 
galaxies of their gas is harassment (JMoore et al.ll 19961 ) . This 
is the heating of a satellite, occurring when a satellite's 
DM sub-halo has frequent high-speed encounters with other 
DM sub-halos, turning low surface brightness discs into 
smaller dwarfs. Although this is usually considered as a 
process occurring in clusters containing thousands of galax- 
ies, ACDM models of the Universe reveal that galaxy sys- 
tems are scaled down versions of clusters and should con- 
tain thousandsofdark matter sub-halos s urrounding galax- 
ies ([Klvpin et al.ll 19991 : iMoore et al.ll 19991 ) . which means ha- 
rassment on smaller scales may be a possibility. 

Furthermore, satellite-satellite mergers could occur 
which would turn l ate-type satellites into early types 
iMcIntosh et al.ll2008h . This however is unlikely to be a sig- 
nificant way of increasing the red fraction, due to the low 
numbers of satellites per host observed. 



gram, we have investigated the radial distribution of satel- 
lites, and how the red fraction of satellites varies as function 
of projected distance, host stellar mass and satellite stellar 
mass. 

We find out of a sample of 3 514 isolated galaxies, 1 426 
which host a total of 2 998 satellite galaxies. Separating the 
galaxies into red and populations, we find 41.2 per cent of 
the satellites are blue with red hosts, 33.3 per cent are red 
with red hosts, 22.2 per cent are blue with blue hosts and 
only 3.3 per cent are red with blue hosts. We find mean 
stellar masses and absolute magnitudes of logj^Q M* — 10.00 
and Mr = -19.96, and log^, A^, = 9.35 and Mr = -19.30, 
for the red and blue populations of satellites respectively. 
The average satellite in this sample is found to be ~ 1/lOth 
of the mass of its host. The mean Alogj^y A^* = 1.18. Our 
main results are as follows: 

(i) Parametrizing the projected density of the satellites 
by E(_R) oc R" , we find best fitting slopes of a ~ —1.0 for 
satellites divided into four host stellar mass bins (Fig. [5]). 
By dividing the satellite hosts into red and blue popula- 
tions we find the satellites of blue, low mass hosts to be 
more centrally concentrated. Power law slopes of q = —0.94 
and a = —1.05 for the satellites of red and blue hosts are 
obtained, respectively. Correcting for interlopers the slopes 
steepen to a = —1.13 and a = —1.38 (Fig. [6]). 

(ii) We find there to be a steady increase in the red frac- 
tion of satellites as projected distance decreases for satellites 
with hosts of stellar mass logjg A^* Js 10.5. The red frac- 
tion of satellites with hosts of stellar masses logj^Q A4 < 10.5 
shows no trend as a function of radius (Fig. [7|. 

(iii) Sub-dividing the satellite sample into host and satel- 
lite mass bins reveals that the host mass has the biggest 
effect on the satellite red fraction. Satellites with more mas- 
sive hosts are more likely to be red. Comparing satellites to 
the small neighbours control sample, we find the red fraction 
is 0.13-0.25 higher for satellites of 11.0 < logj,, A4* sC 11.5 
hosts, 0.06-0.15 higher for 10.5 ^ log^oA^* =^ 11.0 hosts, 
and almost unchanged for lower- mass 10.0 ^ logj^gA^, ^ 
10.5 hosts (Fig.E}. 

The effect of environment, as noted by the difference 
between the satellites and small neighbours, can act over 
large range in high mass haloes and appears to be primarily 
a function of host mass and not satellite mass. These results 
suggest that quenching of star formation by strangulation is 
likely to be the main gas removal process in satellite systems, 
with other tidal effects, having a smaller contributing effect. 
Further study of the radial dependence of the instantaneous 
star formation rates of satellite galaxies using Ha emission 
line measurements, is likely to yield more details about the 
quenching mechanism. 
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